Magnetoresistance of an Anderson insulator of bosons 
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We study the magnetoresistance of two-dimensional bosonic Anderson insulators. We describe 
the change in spatial decay of localized excitations in response to a magnetic field, which is given by 
an interference sum over alternative tunnelling trajectories. The excitations become more localized 
with increasing field (in sharp contrast to generic fermionic excitations which get weakly delocal- 
ized): the localization length £(B) is found to change as £,~ 1 (B) — £ _1 (0) ~ -B 4//5 . The quantum 
interference problem maps onto the classical statistical mechanics of directed polymers in random 
media (DPRM). We explain the observed scaling using a simplified droplet model which incorpo- 
rates the non-trivial DPRM exponents. Our results have implications for a variety of experiments on 
magnetic-field-tuned superconductor-to-insulator transitions observed in disordered films, granular 
superconductors, and Josephson junction arrays, as well as for cold atoms in artificial gauge fields. 
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Transport in Anderson insulators [l|, [|| is crucially de- 
termined by the properties of localized wavefunctions. 
Their structure is very complex, both deep in the insula- 
tor, as well as upon approaching the derealization tran- 
sition, where they develop a multifractal structure Q. 
A particularly important tool in probing the nontriv- 
ial structure of localized states in Anderson insulators is 
magnetoresistance. This is because a magnetic field sen- 
sitively affects the quantum interference which in turn 
influences quantum localization. This effect of the mag- 
netic field has been studied extensively in the past con- 
centrating mostly on non- interacting fermions [J, [H| . 

Recent experiments on disordered superconducting 
films provide evidence for bosonic insulators with local- 
ized electron pairs as carriers @, 0|- These and other 
similar systems feature a giant peak in magnetoresistance 
(MR) 0- 12 1 . This is often interpreted as a crossover from 
bosonic to fermionic transport [131 114| , even though the 
details remain controversial. Bosonic localization prob- 
lems arise also in disordered granular superconductors in 
the insulating regime, in cold bosonic atoms in speckle 
potentials (where artificial gauge fields can mimic a mag- 
netic field) as well as in disordered quantum magnets. 

The predominant mode of transport in disordered in- 
sulators is variable-range hopping of carriers between lo- 
calized excited states [lfj. The spatial decay of wave- 
functions describing these localized excitations deter- 
mines the inelastic hopping rate and thus the resistance. 
At low temperature, the (phonon-assisted) hops become 
significantly longer than the average distance between 
impurity sites hosting the excitations. In this situation, 
one needs to know the wave-function amplitudes at dis- 
tances greater than the Bohr radius of an impurity state. 
At these distances, the amplitude is reinforced by mul- 



tiple scatterings from intermediate impurities whereby 
many alternative paths interfere with each other 0, [H[ . 



A perpendicular magnetic field affects the interference 
of the scattering paths on all length scales and modi- 
fies the localization properties. Interestingly, bosons and 
fermions behave very differently in this respect: while in 
the absence of a field fermion paths typically come with 
amplitudes of arbitrary signs, low energy bosonic am- 
plitudes are positive and thus interfere in a maximally 
constructive way. The magnetic field suppresses this in- 
terference, yielding a strong positive magnetoresistance. 
It exceeds by far a largely opposite effect seen in fermions, 
which arises from a subtle suppression of negative inter- 
ferences 11611 . 



Despite numerous studies of fermionic MR HEMl, 
a full understanding of the effect of magnetic field on 
the large-scale structure of localized wave-functions has 
not been obtained. In this Letter we study the bosonic 
cousin of this problem and show that it is amenable to 
a complete solution. The simplifying circumstance is the 
absence of additional sign-factors in the latter quantum 
interference problem, which allows a mapping to classi- 
cal statistical mechanics of directed polymers in random 
media (DPRM). More generally, our analysis of MR is 
also valid for fermionic problems, provided the interfer- 
ing paths have essentially only positive amplitudes. This 
arises, e.g., in the tunneling through the bottom of the 
conduction band in a solid semiconductor solution [20j |. 
or in fermionic impurity bands with Fermi level very close 
to the band bottom 21| . 



The model - Here we study a model of hard-core bosons 



FIG. 1. The approximation of directed propagation [J maps 
the wavefunction to a directed polymer. The droplet pic- 
ture suggests that traces of localized wavefunctions, or low 
energy polymer configurations, form a string of loops of com- 
peting/interfering paths. Relevant loops of size I have trans- 
verse roughness ~ ^"• =2//3 ). They are rare, being separated by 
a typical distance £ 1+e = £ 2C ~ 2> t- Two competing paths I\2 
are shown, and the loops/droplets they form. 



on a square lattice, 



H = ^(£; - n)c\ci 
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with uniformly distributed on-site disorder in the range 
£( e [— W, W]. We take W = 1 as the energy unit and 
consider weak nearest-neighbor tunneling, t <C W . We 
fix the chemical potential to /z = to study a half-filled 
impurity band. A perpendicular magnetic field is intro- 
duced via the vector potential A = Bx e y , with B being 
the flux per plaquette in units of the flux quantum. 

We now focus on the spatial structure of an excitation 
localized around site i. It is characterized by the residue 
of the pole at to w ti of the retarded Green's function 
Gft{w) = i / °° ' dte iuJt {[ Cj (t),c\(Q)]) [H. Its decay away 
from the site i defines a localization length. Deep in the 
insulating regime, Gj^ can be evaluated using a locator 
expansion [16j |. To leading order in small hopping one 
obtains a sum over all paths T of shortest length Q, 
dist(ij) (cf. Fig. [TJ only right-going steps are allowed) 
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which is closely analogous to the sum over paths for 
fermionic Anderson insulators [l| . In Eq. © each path 
r contributes with an amplitude 
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FIG. 2. Magnetoconductance of fermions and bosons as a 
function of distance iV in a half filled impurity band (fj, — 
0). The linear dependence implies that the magnetic flux B 
changes the localization length £. While it increases slightly 
for fermions, it shrinks much more substantially in bosons. 



by the presence and absence (respectively) of the factor 
sgn(efc) in the amplitudes ©. For bosons, the ampli- 
tudes are all positive for £, = 0. A magnetic field de- 
stroys this complete constructive interference, and thus 
localizes the wavefunction more 3, 18, 2^. In contrast, 
typical fermionic problems [H feature amplitudes which 
vary in sign, depending on the number of sites on the 
path with Si < fi which are occupied in the ground state. 
In this case the dominant effect of a magnetic field lies 
in destroying negative interferences of competing paths, 
which tends to delocalize the wave function slightly. Both 
cases are readily amenable to efficient numerical studies 
via transfer matrices 0, [I3|i which we use below. The 
results shown in Fig. [2] illustrate the opposite trends. 

The relevant quantity for transport is the typical spa- 
tial decay of localized excitations. Therefore one focuses 
on the (typical) magnetoconductance, defined as 0| 



Aa N (B) = cxp (hi[|W(B)/ %( o)|]) , N = dist(zj), (4) 

where the over bar denotes the disorder average. We take 
on opposite corners of a square (cf- Fig. []}. 
The linear variation with distance in Fig. [2] implies that 
at large scales B changes the typical decay rate, i.e., the 
inverse localization length l/£, of the excitations. 

Numerical evaluation - One numerically evaluates 
Sji(B) = S XjtVj (B) (with i as origin) by recursion 



and an accumulated phase §r(B) = J r dr ■ A. On av- S x+ i^ y (B) = V x+ \, v [e 4 * S XiV _i(B) + e l<t>+ S x ^ y+ i(B)](5) 



erage, the larger the excitation energy e^, the faster the 
spatial decay of IS^I [l6|. Henceforth, we focus on low- 
frequency excitations (relevant for transport at low T) 
and hence set lu = E{ = 0. 

Within this "forward-scattering approximation" Q, 
justified for t <^ W, bosons and fermions differ only 



with0± = / r A-(ir, where r± : (x,y±l) — > (x+l,y) are 



straight paths along the lattice links and V x 



Act at (B) evaluated from this varies as B 2 N 3 for small 
(B, N) and shows a sharp crossover to NB 4 / 5 at larger 
fields/distances (cf. Fig. The data for different N is 
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N=200 
N=600 
N=1000 
<i>(x«1) = 0.31 x 1[ 
«(x»1)=0.56x 4 ' 
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It is interesting to note that for weak fields, Eq. J5|) 
admits a continuum limit, where S obeys the equation 



FIG. 3. Scaling of the magnetoconductance, Act, with dis- 
tance N and flux per plaquette, B. The crossover from 
the perturbative regime \\nAa N (B)\ ~ B 2 N 3 to the non- 
perturbative regime |ln Acrjv(B)| ~ NB 4 ^ 5 occurs at N ~ Ib, 
where many successive interfering loops start contributing. 
Inset: change of inverse localization length for N = 200, and 
best fit to the leading two terms in Eq. |QTJ, £ _1 (-B)-£ _1 (0) = 
Cl B 4/5 +c 2 B. 



found to collapse onto a scaling function 

| In Actat (_B) | = AT-V 3 $ (iVB 3 / 5 ) , (6) 

$(x < 1) = bxx w > 3 ■ $(x » 1) = b 2 x<\ 

with b\ w 0.31, £>2 ~ 0.56. This scaling is expected theo- 
retically from the physics of directed polymers (DPRM) , 
as we explain below. 

Mapping to directed polymers - By virtue of the posi- 
tive path amplitudes Sji (B = 0) can be interpreted as the 
partition sum of a DPRM in 1+1 dimensions [25|,[26| with 
random onsite energies In |ej| (at temperature T = 1) and 
ends fixed at sites i and j. Each polymer configuration 
corresponds to a directed path T of the expansion (J5J. 

In low dimensions, DPRM exhibit a pinned phase at 
large scales, as the random potential is relevant under 
renormalization [13, HH- Beyond a characteristic pin- 
ning scale L c (of the order of the lattice scale here), 
the random potential competes strongly with the poly- 
mer's entropic elasticity and induces roughness exceeding 
that of random walks: On longitudinal scales £, typical 
transverse excursions of configurations grow as £^ with 
C > 1/2. A low energy excitation that differs from domi- 
nant configurations on scale £, has typical excitation en- 
ergy E(i) ~ £ 6 , with energy exponent 9 = 2£ — 1 [2{|. In 
1+1 dimensions (MR in 2d), the value C = 2/3 is known 
exactly while C3d ~ 0.62 is known numerically |3lj |. 

When B ^ 0, the polymer configurations acquire com- 
plex weights. Studies of £ and 8 exponents of complex 
DPRM [32| suggest that the scalings of the pinned phase 
do not change with complex weights. In fact, for fermions 
at B = 0, where negative weights are abundant, there is 
numerical evidence that the wavefunctions are still gov- 
erned by DPRM exponents [33M35| . We thus assume that 



D x S = D 2 y S 



V{x,y)S, (7) 

with a (5-correlated random potential term V(x, y) and 
B a =(x,y) = d a — iA a (x,y) being the gauge-covariant 
derivative (in Landau gauge A y = 0). This generalizes 
the Kardar-Parisi-Zhang (KPZ) equation j3(| to the pres- 
ence of complex potentials V — > V + iA x , and may render 
bosonic MR amenable to a field theoretic analysis sim- 
ilar to Refs. [37, 38 1 . However, a rigorous study of this 
modified KPZ equation is not attempted here. 

In DPRM language, the magnetoconductance can be 
cast as a thermodynamic average of the phase factors 
e i*r(B) over polymer configurations, and the ratio of am- 
plitudes Sji takes the manifestly gauge-invariant form: 

~ E r> COS {BA TT , 
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(8) 



Here Er = Sfegr\i m l efe l ^ s * ne cncr gy of configuration 
r, and Arr' is the oriented area enclosed by T and V. 

MR in weak fields - For weak fields or short distances 
one can evaluate Ao~n(B) perturbatively in B. Typical 
loops of linear extent I enclose a flux ~ B£ 1+< > . Of the 
N/l possible independent loops only a fraction ~ £~ e 
interfere significantly, cf. Fig [I] and are thus sensibly af- 
fected by B. As long as N <C £b = B~X+? the dominant 
contribution to Eqn. [8] comes from the largest loops of 
length £ ~JV, which nevertheless enclose only a fraction 
of a flux quantum. This results in the magnetoconduc- 
tance (0} Aa N cx -N- e (BN 1+ () 2 = -B 2 N 3 . Note that 
the roughness exponent drops out of this perturbative 
result. We therefore recover the same scaling as previ- 
ous authors predicted for interfering paths with positive 
weights [3| , even though they assumed random walk scal- 
ing, C = 1/2. However, this coincidence hides the fact 
that typical wavefunctions are less strongly affected by 
B than might be suggested by AaN, since the disorder 
average is dominated by rare events. 

MR in strong fields - For N>£b, DPRM scalings show 
more clearly in the magnctorcsponsc. The dominant con- 
tribution to Ac tv comes from reduced interference in 
loops of length £b, each of which decreases Acta? by O(l). 
Larger loops contribute similarly, but their probability 
to interfere significantly decreases as £~ e . On the other 
hand, smaller loops, albeit more abundant and likely to 
interfere, enclose a small fraction of a flux quantum, and 
thus have a negligible effect. The contribution from loops 
of size £b gives rise to an extensive ln(AcrAr) proportional 
to the density of significantly interfering loops, 



the DPRM exponents hold for finite fields as well. 



lnACTjv _ A (\ 

This is equivalent to a reduction of the inverse localiza- 
tion length by B 4 / 5 in 2d. In 3d the same arguments 
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FIG. 4. Hierarchical droplet model: At each level of the hi- 
erarchy, a parent loop C (composed of a dominant and sub- 
dominant branch) is split into four subloops, two forming the 
dominant branch (£'1,2, thicker line), and two forming the 
subdominant branch (>C" 2 , thinner lines), cf. Eq. (|10[) . The 
parent levels are indicated by dashed lines. The dots indicate 
the splitting into two successive loops at the next level. 



apply, with an exponent 2£/(l + £) ~ 0.765. Both ex- 
ceed the value 2/3 obtained upon neglecting pinnin g an d 
assuming random walk scaling with £ = 1/2 0, EH, Hil . 

So far we have discussed the leading scaling with mag- 
netic field. However, the numerical data show small sub- 
leading corrections (cf. inset of Fig. [3]). Those are in- 
deed to be expected from spatially overlapping loops. 
To understand their effect, we introduce a hierarchical 
model which incorporates the essential ideas of droplet 
theory for directed polymers [H, [2^] . At a given length 
scale L, the polymer has typically a preferred set of con- 
figurations, which compete with alternative, subdomi- 
nant sets of paths. The leading subdominant family of 
paths has a higher free energy by 0(L ) and wanders 
off the dominant configuration by , enclosing a typical 
loop area 0(L 1+Ci ). This pattern repeats at all length 
scales. We simplify this phenomenology by considering 
a model where loops and alternative paths are restricted 
to lengths Lk = N2~ k where N ^> 1 is the fixed dis- 
tance between endpoints. Each parent loop of size Lk is 
composed of a dominant and a subdominant set of paths, 
each being made up of two successive loops of size Lfc+i, 
cf. Fig. |4] We define the propagation amplitude over the 
distance N recursively. For a parent loop C at level fc we 
encapsulate DP scaling by defining the amplitude 



S 



-fcLijacBLt 



• cfe+1 cfc+1 

O r ,i Of,, , 



(10) 



where C x 2 and CJ[ 2 are the child loops along the dom- 
inant and the subdominant path, resp. fc > and 
ac arc random variables of order O(l), with a proba- 
bility density p(fc, a c), assumed to be i.i.d. for all loops 
C. The recursion is closed by setting all = 1 for fc 
with Lk < (-b [40]. The magnctoresistance is defined as 

Aa N = ln(|S O jv(B)/S(w(0)|). 

This model has elements in common with the hierar- 
chical lattices analyzed in Rcf. 41- However, here we 



explicitly include the known scaling of excitation ener- 
gies and areas of loops. The latter is necessary to discuss 
physically meaningful magnetorcsponse. Note that sig- 
nificant interference between the paths C and C" (as 



given in Eqn. [TOf occurs only for rare 'active loops' £ for 
which f c L e < 1. 

The perturbative scaling Acr/v ~ B 2 N 3 is easy to ob- 
tain in this model [42j] . In the non-perturbative regime 
(N S> £b 3> 1), using that active loops are sparse, one 
can expand In Ao~n{B) in p owers of the density of active 



loops of linear size Ir. [43 



\uA<7n(B) „^_r 1-.S- ,-,-22- 



L (ii) 



where the constants a depend only on the distribution 
p(/, a) . Similar to the Mayer cluster expansion for a 
system of interacting particles, one obtains a term of 
O (£? (1+ " £ ' ) A 1 +c)^ by collecting contributions with exactly 
n active loops. The leading coefficient c\ is positive defi- 
nite, and we found ci > 0, independently of our choice of 
the distribution p(f, a). Subleading terms due to interfer- 
ing loops thus enhance the negative MR of bosons. This 
may explain a similar effect seen in the numerical data 
on the original lattice (inset of Fig. [3]), where a fit yields 
Ci ss 0.34, C2 ~ 0.67. Hence, In Act appears to follow a 
power law with slightly larger exponent than 4/5. 

Experimental consequences - In transport through 
variable-range hopping, at fixed T, the resistance de- 
pends on the localization length as R(£) = pexp(A/£ Q ), 
with a = 1 /i (with Coulomb gap) and 2 /3 without (Mott's 
law in d = 2) [l5|. According to © a perpendicular mag- 
netic field reduces the bosons' localization length £ as 

l/£(B)«l/£ + A[l/a(B), B min <B. (12) 

where B min = [^ln(i?(0)/p)] _5/ ' 3 is needed for l B to be 
shorter than the typical hopping distance. To lowest or- 
der this effect increases the resistance by the factor 



R(B)/R(0) = [R(0)/p] 



a£A[l/f] 



Bmin < B. (13) 



For B <C 1, the exponent is aci£B 4 / 5 . For B -> 1, 
it receives subleading enhancements, reaching values as 
big as 0.3a£ Q, cf. Fig. [3] As resistances up to 
R(0)/p ~ 10 6 are measurable, and localization lengths 
£ < 2 are expected to be within the regime of appli- 
cability of forward scattering (loops being sufficiently 
suppressed) our theory predicts strongly positive MR of 
bosons, with enhancement factors of up to two orders 
of magnitude, within the theoretical and experimental 
limits. These effects are even stronger when resonances 
are suppressed 44|. A further enhancement of MR is 
expected in the critical regime where loops must be in- 
cluded. In contrast, the analogous fermionic problem 
exhibits negative MR, which moreover reaches a much 
smaller maximal amplitude, cf. Fig. [2] The importance 
of the bosonic MR makes it likely to be a key ingredient 
in the MR peak observed in superconducting films with 
preformed pairs 0. Finally, it would be interesting to 
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probe for the predicted magntcorcsponse and its sensi- 
tivity on quantum statistics using cold atoms subjected 
to artifical gauge fields. 
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SUPPLEMENTARY MATERIAL FOR BOSONIC 
ANDERSON INSULATORS IN A MAGNETIC 
FIELD 

I. HIERARCHICAL LOOP MODEL 

Here we analyze in more detail the hierarchical loop 
model defined in the main text. This model imple- 
ments the ideas of the droplet picture in an analytically 
tractable and mathematically precise way. Our focus will 
be on the analytical calculation of magnctoconductancc 
in the perturbative and the non-perturbativc regimes. 
However, we have also studied the crossover between the 
two regimes numerically. We found the crossover, the 
asymptotic power laws and subleading corrections to be 
very similar to those observed in the full lattice model of 
forward-directed paths. This suggests that the hierarchi- 
cal droplet model captures indeed most of the relevant 
physical ingredients of magnctoconductance. 



A. Models 

Imposing the scaling of individual droplet degrees of 
freedom actually does not fully specify a hierarchical 
droplet description, but leaves some freedom in the def- 
inition of the model. The resulting models differ in the 
way they treat correlations between energies of spatially 
overlapping droplets. As we will see this translates pri- 
marily into differences in the numerical coefficients of 
subleading terms. 



1. Normalized recursion 

We first discuss a different version of the hierarchical 
construction from the one in the main text. We define 
it by iterating the following recursive construction from 
the largest scale N down to the lattice scale, the loops or 
branch segments having lengths Lk = 2~ k N for < k < 
K = [log 



Sc. 



= 1, if 
S C[ S 



Sc = 



W C {B) = exp 



exp 



L c = 2- K N, 
'c- 2 +Wc(B)Sa { S ai 



1 + 
-fL'c 

-fL e c 



W C (P) 

+ iBa c L c +C 



(14) 
(15) 



ia c 



Lc 

£b 



i+C 



which differs by the normalizing factor 1 + Wc(0) from 
Eq. (10) in the main text. We have defined £b = 
.B _1 /(i+C) and have dropped the explicit dependence of 
Sc on B. Note that the normalization factor in the de- 
nominator in Eq. (|14|) ensures that Sc(B = 0) = 1. 
Therefore fL c is precisely the free energy difference 



between the leading and subleading branches of paths, 
which this model treats as independent from loop ener- 
gies at smaller scales. In weak fields the magnetic field 
response will be insensitive to the precise value of the 
small scale cutoff, L m [ n = N2~ K , as long as it is much 
smaller than the relevant magnetic length, L nlln <C ^b- 
Indeed, up to small corrections, Sc{B) as 1 for all loops 
with L c < £b- 



2. Non-normalized recursion 

The above model assumes that free energy differences 
between a dominant and subdominant branch are in- 
dependent of the energies (and thus interferences) on 
smaller scales along those branches. A more realistic 
model should take into account that if positive interfer- 
ences occurred along a branch, the resulting " free energy" 
of the branch is statistically smaller than if the interfer- 
ences were negligible. Such effects can be built into a 
hierarchical construction by modifying the recursion to 



Sc 



1, if L c 



Sc = S C[ S C , + W £ (B)S C »S C >>, (16) 

with the same weight factor Wc(B) (fl5)) . but dropping 
the normalization. In this case £x(0) is not normalized to 
1 at all length scales. Instead, the explicit contribution to 
free energy difference between two branches, fL^, is now 
supplemented by an extra contribution coming from the 
sum over paths at smaller scales. This non-normalized 
recursion follows a similar hierarchical construction by 
Derrida and Griffith [4l|. Those authors assigned to 
each loop random energies or signs, that however did not 
scale with the level of the hierarchy. This generated ran- 
domly fluctuating free energies, with a free-energy expo- 
nent which is only slightly smaller than the value 9 = 1/3. 
In this version of the recursion, we retain the spirit of the 
Derrida-Griffiths approach with the difference that we in- 
troduce the DP-scaling by hand through the free-energy 
fL c , and do not consider random signs in the recursion 
relation. 



B. Magnetoconductance 

We now study the magnctoconductancc of the above 
models, 



In l\(Tn(B) = In \S C={Q N}(B)\ - In |S £={OJV} (0)|j(17) 



where [• • • ] denotes the average over the set of reduced 
free energy and area variables, {h = (/, a)}. We assume 
the two variables associated with each loop to be inde- 
pendent and identically distributed. 



P({h}) = Y[p(fc,ac)df c da c . 



(18) 



7 



The product runs over all loops C, the support of p being 
{/, a} € [0, oo ) x (—oo, oo). However, as we will see, only 
the values of p(f = 0, a) = p a (a) will enter the analytical 
results. For quantitative calculations, we will assume a 
simple Gaussian form, 

p(f = 0,a) = Pa (a) = p eXp[ Zf /2a2 ° ] . (19) 

V Z7T C(o 

The (non- normalized) density p a {a) is a free input pa- 
rameter of the hierarchical models. More realistic densi- 
ties could be determined by studying the distributions of 
loop areas in the full lattice model. 

Let us now analyze the magnetoconductance, 

F({h}) = [ln\S 0N (B)\-]n\S 0N (Q)\]({h}) (20) 

as a functional of the disorder realization ({h}). F can 
be viewed as the free energy difference between a directed 



polymer with i?-induced complex weights and one in zero 
field, where all weights are positive. 

In typical disorder realizations most loops do not play 
a significant role in modifying the interference of alterna- 
tive tunneling paths. A loop L is involved significantly 
only if fc < L^ 6 , in which case we refer to it as 'active'. 
Large active loops are dilute, while small ones are more 
abundant, but contribute very little to magnetoconduc- 
tance. One can thus expand F in the spirit of a droplet 
or virial expansion into a sum of terms Vk , which involve 
an increasing number k of spatially overlapping loops, 

F({h}) = V1+V2 + V3 + ... 

= E £ F c (h Cl ,...,h Ck ) . (21) 
fe>l {Cijt.. .^c k } 

The sums are over all (non-ordered) sets of distinct loops. 
The decomposition in Eq. (|2Tj) is exact, given that the 
connected functions F c are defined recursively as 



F c (h Cl ) = F ({h c \f c ^ Cl ^00}), 
F c (h Cl ,h C2 ) = F ({h c \f c ^,c 2 -> 00}) - F c (h Cl ) - F c {h C2 ), 



(22) 
(23) 



fc-i 



F c (h Cl ,...,h Ck )=F({hc\fc^ooV£t{C 1 ,...,£k}})- ^ ]T F c (h c , , h a J. (24) 

m=l {C' 1 ^...^C' m }c{Ci,...,Ck} 



The subtraction of the disconnected terms in Eqs. (|23l24p 
ensures that F c tends to as one of its free energy ar- 
guments becomes large, ft — > 00, which turns the corre- 
sponding loop inactive. It is also easy to verify that F c 
vanishes, unless the loops associated with its arguments 
belong to a single spatially entangled cluster. This fol- 
lows immediately form the fact that disconnected sets 
of loops contribute additively to In \Sqn(B)\. This clus- 
tering property ensures an extensive result in the large 
distance limit, N ^> £b (for every order of the expansion 
Vk ~ N), i.e., we must have 

In \S 0N (B)\ - In |5 jv(0)| = -A^ 1 )^ + o(N), (25) 

where the coefficient A(£ _1 ) is expected to be self- 
averaging. As the notation suggests, this coefficient rep- 
resents a correction to the inverse localization length 

The disorder average is carried out term by term. 
Thereby, the disorder variables, especially fc, take the 
role of relative positions of particles played in the clus- 
ter expansion of gases. The role of a low gas density is 
played by the small likelihood of large loops to be active. 
The term Vk in the expansion (I21[) captures the inter- 
ference contribution from exactly k active loops, similar 



to droplet expansions at low T in related disordered sys- 
tems. [45j-|47| This is akin to the virial expansion, which 



corrects the ideal gas behavior by summing ?i-particle 
contributions at order n k in an expansion in the density 
n. 

The various contributions to Vk can easily be repre- 
sented graphically by enumerating all spatially connected 
sets of k loops, and summing over their sizes, see Figs. [5] 
and [6] . 



C. Evaluation of leading terms 

1. 1 st order term 
The first term in Eq. (|2ip can be rewritten as 



V 1= £ln 



l + W c (B) 



l + Wc(0) 



(26) 



2 ^ 

c 



In 



4e-^ L ° sin 2 ( ZfBL c +C 



(l + e -f* L i 
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for both the normalized and the non-normalized recursive 
definitions of the model. 

Reorganizing this as a sum over looplcngths £k = 
N2~ k , and performing the disorder average, we find 

-aj 



K 



1 - 



i+C 



cosh' 



(27) 



2. 2 nd order term 

The second term in the droplet expansion, V 2 , picks up 
contributions from disorder realizations where two active 
loops spatially overlap. This can occur in two distinct 
ways, c.f., Fig. [5] either (I) the smaller loop is part of 
the dominant; or (II) part of the subdominant branch of 
the larger loop. Let us refer to the bigger and smaller 
loop as C\ and C 2 , respectively, with lengths L\^ 2 . 

The following expressions apply to the normalized 
model. The discussion of differences for the non- 
normalized version will be discussed further below when 
we evaluate the terms. Denoting Wi{B) = Wci(B), V2 
can be written as 



V 2 



£1X2 

Li>L 2 



[v 2 (J) (A,£ 2 ) 



V [ 2 I] {Ci,C 2 ) 



(28) 



where 



Vf ) (£i,£ 2 )=ln 



1 + W 2 (B) 
1 + W 2 (0) 



W^B) 



1 + wMo) 



(29) 



In 



l + Wi(B) 



V ( 2 II) (£ 1 ,£ 2 ) = \n 



1 + Wi(0) 

, 1+W 2 (B) 
r 1+W 2 {0) 



-In 



W 2 (B) 



1 + W 2 (0) 



W^B) 



In 



1 + Wi(0) 
l + Wi(B) 



1 + Wi(0) 



(30) 



Taking the disorder average and writing (|28|) as a sum 
over loop lengths we have 



K K 



N 



1 



W 1 (B)(Z 2 + Z 2 7-2) 



(l + W!(B)y 



ki—0 k2>ki 

where Z 2 = (1 + W 2 {B))/{1 + W 2 (0)) 



a i,2,/l,: 



(31) 





FIG. 5. Graphic representation of the first two virial terms in 
Eg. 1211 Left: Vi is a sum over all loops C, composed of a dom- 
inant (thick line) and subdominant (thin line) branch. Right: 
The two contributions to V2 arise form spatially overlapping 
loops of length Li > Li. The two cases distinguish whether 
the smaller loop belongs to the dominant (I) or subdominant 
(II) branch. 

the other hand, we will see that small scales contribute 
negligibly to magnetoconductance as long as B <C 1, so 
we can concentrate on Lc > is- Thus, for each vari- 
able / in the disorder average, we can safely approximate 
/ dfdap(f,a)... w / dfdap a (a)..., cf. Eq. dHJ). 

Introducing Fi = fd L C . , the disorder-average of the 
n-th virial term becomes 



{£!#...#£„} V=l 



In, 



(32) 



= 1 



J « = II [J dF i J Pa(ai)d ai j F C (H U H 2 , ■ ■ ■ ,H n ), 



using the notation Hi = (Fi , 



II. SCALINGS IN THE DROPLET EXPANSION 



3. Higher order terms 

For sufficiently large loops, Lc 3> 1, the disorder aver- 
age simplifies. Indeed only very small values of fc { are 
relevant, since the connected functions F c fall off rapidly 
when one if its arguments fc i is larger than L^ 6 . On 



A. Weak fields: j5iV 1+c < 1 

For weak fields one can expand /„ in the enclosed 
fluxes, the result being dominated by the largest scale 
N. Expanding Eq. (|2"T|) in B, and integrating over the 
rcscaled F variables, we find the leading contribution to 
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FIG. 6. Graphic representation of the third order terms V3 in 
Eq. dHJl with Li (blue) > L 2 (red) > L 3 (green). 

the magnetoconductance lnAcrjv(B) 



1 r K 

Vi w-j / p a (a)a 2 da 5 2 iV 

k=0 

f p a {a)a 2 daB 2 N 3 , 



(33) 



which is negative, as expected for bosonic magnetocon- 
ductance. Likewise, one can check from (|3Tj) . that V2 ~ 
0(B 2 N 2 ( 1+ Q- 2e ). More generally one finds that higher- 
order terms are suppressed by the prefactors Yl2=i ^1° 
in Eq. (f3"2"f with Li ~ AT, which leads to the subdominant 
scaling V fc - 0(B 2 AT 2 ( 1+ f)- fce ). 

Note that the leading scaling (|33|) is independent of the 
wandering exponent, by virtue of the relation 6 = 2Q — 1. 
One therefore obtains the same scaling as in a non- 
disordered case, for which the exponents ( = 1/2,9 = 
hold. However, we stress that in the disordered case the 
result (f3"5")) arises as a result of disorder averaging, which 
masks some of the physics. The distribution of Vi is wide, 
and the average (|33|) is dominated by a few rare dis- 
order configurations. The latter occur with probability 
~ N~ e , but contribute a large Vi ~ B 2 N 2{l+ ^\ while in 
most other realizations the wavefunctions are much less 
affected by quantum interference. (48j 



B. Strong fields: BN 1+( > 1 

The perturbative expansion holds only for weak fields 
for which the distance between end points is smaller than 
the 'magnetic length', TV <C £b- For stronger fields, the 
dominant contribution comes from loops at the scale £b- 
To see this, let us approximate the sum over discrete loop 
sizes in (|27|) as an integral, J2k ~ V m (2) / di/£, 



Vi 



N 



21n(2) 



di 



In 



1 - 



sin 2 (f BP-+<) 



cosh" 



■a,f 



•(34) 



Rescaling the free energies and changing variables to u : 
£/£b, we obtain 



Vi 



-ci 



N 
JT+e 



- Cl NB^< 



-ci 



NB 4/5 , (35) 



with the numerical coefficient 

1 f°° du 
Cl ~ 21n(2)./ u 2 +<> 

dap a (a) In 



dF x 



2 {au 1+ </2 



cosh 2 (F/2) 



(36) 



For the particular choice (|19|) for p a {d) (with po = ao = 
1) we find 



ci w 0.86. 



(37) 



Note that the dominant contribution comes indeed from 
u = O(l), i.e., from loops of size Ib- We have extended 
the limits of the M-integral to and 00, as it converges 
rapidly on both sides. This result implies a leading cor- 
rection to the localization length as 

AOT 1 ) = ci^ttc. (38) 



C. Subleading corrections 

In the non-perturbative regime, subleading corrections 
are interesting to analyze in more detail, as they correct 
the leading behavior (|38p . As we shall see below, there 
is a direct correlation between the order of a term in the 
virial expansion and its scaling with B, which justifies 

using the virial expansion in the first place. Roughly 

s 

speaking, each loop contributes a scaling factor of B x +( 
on disorder averaging, causing the n-th virial term with 
n spatially overlapping loops to contain as many such 
factors. 

A more precise formulation follows. We begin with 
the second order term V2- As before, we can rewrite V2 
as a sum over pairs of loop lengths, Lj = 2~ J A r and 
Lh = 2~ k N, and take the continuum limit of the discrete 
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sums over j and k 



N 



3,k 
k>j 



N 



N d£ 1 f^dt 



(39) 



the appropriate probability distributions 

dfiPf(h) / dhp f (h) 
Jo 

da 1 p a (a 1 ) / da 2 p a (a 2 ) 



(ln2) 2 J 1 l x J x I 



Substituting u\ = and u 2 = j^, we obtain the 
sublcading correction 



V 2 = -c 2 - 



N 



(1+29) 



Here, (£i,i 2 ) and V^ 5 (4,40 are given by 

Eqs. (|29|30jl (loop d being referred to by its length li). 
The overbar denotes the average over /i j2 and ai j2 with with the numerical coefficient 



4<:-i 

-c 2 NB^+T. 



(40) 



B 



C 2 



(In 2) 



dui 



o Ml 



o u 



2+6 



daip a (ai) / da 2 p a (a 2 ) 



x In 



e -Fi+mitij 



^ _|_ e -F 2 +ia 2 U2 + 



^ + e -F 1 +iaiu 1 i 



1 + e 



-F 2 



(41) 



Note that the integrals converge both for u\ i2 — > and 
1*1,2 — > CO- 

We have computed the (i*i, -F 2 , &i, a 2 )-integral in 
Eq. ([4Tj) using Monte-Carlo sampling, followed by nu- 
merically carrying out the (ui, u 2 )-intcgration. With the 
density p a given in (fT^|) , c 2 turns out to be negative. 

To the first subleading order we find the correction to 
the inverse localization length as 

A(r') = -ciB 4 /' 5 ^1 + J^5 1/5 + 0(B 2 / 5 fj . (42) 

Note that the subleading corrections vary slowly with B 
and thus are expected to affect fits of the magnetocon- 
ductance to a simple power law B 1 . Indeed defining an 
'effective exponent' as 



7 



dh\B 



1 



5d/ca + SVs 



0(£ 2/5 )(43) 



one expects to see apparent exponents that deviate from 
the asymptotically exact value 4/5 for any small but finite 
B <C 1. The sign of the correction depends on the relative 
sign of c\ and c 2 . 

The numerical data obtained for the full lattice model 
is consistent with a positive correction to the exponent, 
cf. inset of Fig. 3 in main text.. However, the normalized 
hierarchical model predicts the opposite sign. We believe 
that this qualitative difference is due to the fact that the 
normalized recursion neglects correlations of free energy 
differences at different scales, as explained above. 

A more realistic model, which builds in such correla- 
tions was given in Eq. (|16p . where the normalizing factors 
are dropped in the recursive definition of path weights. 
The expression for V 2 is easy to derive in this case as 
well, 



c 2 



V. 



In 



dui 



(ln2) 2 7 A uY 



du 2 

..2+6 



dF, 



]_ _|_ e -F 1 +ia 1 u\ + ^ _|_ e -F2+M2«2 +<: 



if" = In 



,-F, 



o-F 2 



1 + e -F 1+ ioi«; +c (1 + 



3 -F l(1 



-F 2 ) 



dFo 



daip a (at) / da 2 p a (a 2 ) 



In 



1 — 1 1 14 
I _|_ g-Fi+iaiu.^ 



-In 



,-F, 



1 + e 

2 _|_ g-Fi+iaiuJ 



1 



V, 



(I) 



In 



e ~F 2 +ia 2 ul +i 



1 + e 



-F 2 



V. 



(44) 



I 

where the lower cutoff A is a number < 1, ensuring that the recursion ends at Lk = A£b- Wc used A = 0.125 for 
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the numerical computation of c 2 below. This cutoff is 
required since the ^-integral in Eq. (|44[) does not con- 
verge at infinitcsimally small length-scales. This reflects 
the fact that the Sc (and thus the loop free energies) 
have a non-trivial distribution already in B = 0, due to 
interferences at small scales L < I5. This distribution 
cannot be captured easily by the virial expansion. In- 
stead we have to introduce a small scale cut-off at some 
fixed length scale Lk < Ib- We can safely assume that 
the small scale interference is incorporated into the free 
energy differences at that smallest scale. Thereby we 
rely on the fact that smaller loops enclose negligible flux 
and thus do not contribute significantly to magnetocon- 
ductance, nor affect much the free energy distribution at 
small scales. Finally, this prescription leads to a similar 
virial expansion in powers of 5 1 / 5 , however with different 
coefficients Ck>i- 

The integral flU} yields c 2 «4.9xl(T 2 . This has the 
same sign as ciand thus leads to an "effective exponent" 
which is bigger than 4/5, and thus comes closer to the 
phenomenology observed in the full lattice model, as one 
may expect. As mentioned before, the various definitions 
of the hierarchical construction only affect the coefficients 
of the subleading terms in the virial expansion. 

D. Effect of small denominators and resonances 

The quantitative effect of the subleading terms is of 
course non-universal, as are the coefficients C\ y2 . A varia- 
tion of such effects is actually also found in the full lattice 
sum of forward-directed paths. It may seem dangerous 
to evaluate path sums of products of denominators which 
can become arbitrarily small. While the logarithmic av- 
erage of such sums is mathematically well-defined, it is 
known that backscattering and self-energy effects, or a 
Coulomb gap in the density of states, reduce the influ- 
ence of such resonances. For this reason the toy models 
considered in the earlier literature 0, EH have restricted 
themselves to finite denominators. 

Numerically evaluating the sum over all paths as given 
in the maintext, without restricting the occurrence of 
resonant denominators, we found effective exponent of 
the order of 7 « 0.88. However, the deviation from 4/5 
turned out to be much smaller for a toy model where 
we restricted onsite energies to the interval [1/2, 1]. It is 
thus suggestive to attribute the stronger deviations with 
resonances included to an enhanced value of c 2 . 

E. Higher terms in the droplet expansion 

It is not difficult to write down the disorder-average of 
the higher order terms in Eq.[2T]as appropriate integrals. 



One can check that the generic term Vk varies as 

= c k NB^r (45) 



To illustrate the procedure, we give the diagrams con- 
tributing to V3 in Fig. [6] The corresponding expressions 
for the connected terms are given below. Subscripts 1, 2 
and 3 denote three loops with lengths L\ > L2 > L 3 . For 
brevity, we only consider for the normalized model and 
give the connected terms in V3 as V3 {B) — V3 (B = 
0), where 



V^\B) 


= In 


1 - 


- Wi - 


VW 2 + 


W 3 + W 2 W 3 \, 


V^\B) 


= In 


1 - 


-Wi - 


VW 2 + 


w 3 \, 




= In 


1 - 


-Wi - 


h!F 2 + 


w 2 w 3 \, 


vr\B) 


= In 


1 - 


- Wi - 


YW 2 + 


WiW 3 \, 


VP(B) 


= In 


1 - 


-Wi - 


hWiW; 


i + WiW 3 + WiW 2 W 3 


vr\B) 


= In 


1 - 


-Wi - 


hWiW; 


i + WiW 3 \, 


v r i] {B) 


= In 


1 - 


-Wx - 




i + WxW 2 Wa\, 



where W\ = Wd{B). Continuing along these lines, 
In Act at (B) can be calculated to any desired order at a 
given field B. 



F. Remarks on fermions 

It might be interesting to generalize the hierarchical 
model to the case of fermions. Since the locator expan- 
sion yields path amplitudes with positive and negative 
signs, it would seem natural to include random signs sc 
in a hierarchical droplet model. However, several sub- 
tleties may need further modifications to capture the de- 
tails of fermionic magnetoconductance. For example, a 
weak field can have a significant effect on small loops 
whose branches have nearly opposite amplitudes. This 
may reflect in a non-trivial dependence of free energy 
costs fc on B, which may enhance subleading corrections 
and potentially even change their exponent. It is possible 
that the observed effective fermionic exponents 7 < 4/5 
in the non-perturbative regime are due to such effects. 
More detailed investigations are necessary to clarify these 
issues. 



